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Abstract. In this paper we prove the existence of smooth solutions to fully nonlinear and non- 
local parabolic equations with critical index. The proof relies on the apriori Holder estimate for 
". advection fractional-diffusion equation established by Silvestre 1. 11 J . 

o 

(N 

^ ■ 1 . Introduction and main result 

In this paper we are interested in solving the following fully nonlinear and nonlocal parabolic 
equation: 



(N 



P^'. dfU = F{t,x,u,Vu,-(-A)"u), u(0) = (fi, a e (0,2), (1.1) 

; where F(t, x,u,w,q) : [0, 1] x R'' x R x R'' x R ^ R is a measurable function, and (-A)? is the 
^ ■ usual fractional Laplacian defined by 

(-Ayiu = r-\\-rru), ms^cr^x 

where !F denotes the Fourier's transform, S(W') is the Schwartz class of smooth real-valued 
^ ! rapidly decreasing functions. 

> I Recently, in the sense of viscosity solutions, fully nonlinear and nonlocal elliptic and par- 



abolic equations have been extensively studied (cf. El [TOl [31 13, etc.). In [|4J, Caffarelli and 



00 ■ Silvestre studied the following type of nonlocal equation: 

lau(x) := sup inf I c'^' + b'^ ■ Vu(x) + [ [u{x + y)- uix)]a'j (y)\yr''-" dy] = 0, 

where a 6 (0, 2), /, j ranges in arbitrary sets, c'^ e R and b'^ e R'^, the kernel a'\y) satisfies 

a'^(y) = a'\-y), Uq < a'\y) < Ui. 

^ ' This type of equation appears in the stochastic control problems. In p|, the extremal Pucci 
■ operators are used to characterize the ellipticity, and the ABP estimate, Hamack inequality and 
interior C^'^-regularity were obtained. In IfTTTl . Silvestre studied the following nonlocal parabolic 
equation with critical index a = \: 

dtU = I\U, m(0) = ip, 

and established C^'^-regularity of viscosity solutions. In particular, the following first order 
Hamilton-Jacobi equation is covered by the above equation when H is Lipschitz continuous: 

d,u + H{Vu) + {-A)^u = 0. 

In [|9]|, Lara and Davila extended Silvestre's result to the more general case, and in particular, 
focused on the uniformity of regularity as a ^ 2. 

However, it is not known how to solve the fully nonlinear and nonlocal equation (11.11) in 
Sobolev spaces. Let us fix the main idea of the present paper for solving (|l.ll) . Assume that F 
does not depend on u. Taking the gradient with respect to x for equation (11.11) . we have 

dtVu = -(5„F)(-A)?Vm + iV„F)V^u + V^F. 
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We make the following observation: 

-i-Ayiu = (-A)^divVM = 11" ■ Vu, 

where 'R" := (-A)~div is a bounded linear operator from Bessel potential space to 
provided > 1. If we set w := Vm, then w satisfies the following quasi-linear parabolic system: 

drw = -{d,jF){w, K^wX-AY^w + {V,,F){w, n"w)Vw + (V.,F)(w, n"w). (1.2) 

It is noticed that the classical quasi-geostrophic equation takes the same form (cf. [I6l[5l[8]|, etc.): 

5,0 + (-A)t0 + !Re. V0 = O, !R:=V^(-A)-i 

Assume now that one can solve equation (11.21 ). then it is natural to define 



u{t,x) := ip{x) -\- I F{s,x,w{s,x),'R"w{s,x))(ls. 
Jo 

Thus, if one can show 

Vm = w, (1.3) 

then it follows that 

c 

u{t, x) = (fix) + I F(s, X, Vu(s, x), -{-A)^u(s, x))ds. 
Jo 

For solving equation (II. 2L we shall use Silvestre's Holder estimate [fTTTl about the following 
linear parabolic equation: 

dfU = -a(-A)^u + b -Vu + f. (1.4) 

For proving (11.31) . we need to solve a linear equation like 

dtU = a(-A)'^nu + b-iVu-(Vuy), (1.5) 

where □ := divV - Vdiv is a symmetric operator on L^iW^; W') and 

{ou,u)2 = -||Vm||2 + ||divM||2. 

Notice that in one dimensional case, = 0. 

In this work, we mainly concentrate on the critical case a = I and prove the following result: 

Theorem 1.1. Assume that dqF > ao > and for some kq > 0, 

\F(t,x,u,OM<Ko{\u\ + iy, (1.6) 

and for any R > 0, 

F 6 L'»([0, U-C^iR' xB'^xB^x 5^)), (1.7) 

d,F, V,,F 6 L~([0, 1]; Ci(R' xB'^xBix R)), (1.8) 

5„FeL'"([0, l]x]R^x5j;XR^xR), (1.9) 

where 5^ denotes the open ball in R^ with radius R and center 0; and for any j e N and 
R > 0, there exist Crj > 0, Yrj > and hj{j e {O n L°°)(R'') such that for all (t, x, u, w, q) e 
[0, l]xR''x5^xR^'xR, 

\V{F{t,x,u,w,q)\ < C«j|w|(|wr + D + V/^), (1-10) 

where jr^x = 0. Then for any initial value ip 6 U°° := n^. ^U'^''', where U'^''' is defined by l\2.4\) 
below, there exists a unique u e C([0, 1]; U'") solving equation ALU with a = \. Moreover, 

sup \\u(t)\\co < £""(11^1100 + Ko). 
fe[0,l] 



Remark 1.2. Let A{q) e C°°(R) have bounded derivatives of first and second orders and ^^A 
be bounded below by ao > 0. Let H e C'^iW') and f e ^(IR) satisfy \f{u)\ < kqQuI + 1). Then 

F{t, X, u, w,q) := A(q) + H(w) + f{u) 

satisfies all the conditions ^i.(5D - 47770l) . 

In the subcritical case a e (1,2), when we adopt the same argument described above to solve 
the fully nonlinear equation (11.11) . there are two difficulties occurring: on one hand, we need to 
prove a stronger apriori Holder estimate for equation (11.41) (see Theorem 12 .41 below) 

\u(t,x) - u(t,y)\ 
sup sup ■ < C, 3/3 e (a -1, 1), 

fe[0,l] .v^y \x-yr 

where C only depends on the bounds of a, b, f and m(0); on the other hand, for a e (1,2), it is 
not known whether the uniqueness holds for equation (11.51) in the class of smooth solutions. In 
the case of a e (0,1], this problem can be solved by observing divDa = (see Lemma |5TT] ). 

In the supercritical case a e (0, 1), it is well-known that there exists an explosion solution for 
one-dimensional fractal Burger's equation (cf. EITTl). Nevertheless, from the proof of Theo- 
rem [TTTl one can see that the approach also works for the following fully nonlinear equation: 

a 

dfU = F{t, X, u, -(-A) 2 u), m(0) = if. 

The paper is organized as follows: In Section 2, we prepare some notations and recall some 
well-known facts for later use. In Section 3, we solve the linear equation in Sobolev spaces. In 
Section 4, we prove the existence of smooth solutions for the quasi-linear nonlocal parabolic 
system. In Section 5, we give the proof of Theorem [TTTl 

2. Preliminaries 

Let No := N U {0}. For p e (1, oo) and e No, let W^'' be the completion of S{W^) with 
respect to the norm 

ll/IUp:=Xill^'/IU 

k=0 

where V*" denotes the fc-order gradient; and for < /3 integer, let W^'^' be the completion of 
S(R'') with respect to the norm 

where for a number /3 > 0,[J3] denotes the integer part of /3 and {/3} = /3 - [J3]. It is well-known 
that for k^m,0€ (0, 1) with il-e)k + me^M (cf. [JA p. 185, (2)]), 

(W'^'P ,W"''P)e,p = w^i-^^^^'"^-'', (2.2) 

where (■, Oe.p stands for the real interpolation space. For t e [0, 1], write Y^'' := L^([0, W''''') 
with the norm ^ 

and let xf'' be the completion of all functions u e ^"([0, t]; S{W^)) with respect to the norm 

\\u\\^k.p := sup ||M(i')||/:-l,p -I- \\u\\.^k.p + \\dtU\\^k-Up. 
' se[0,t] ' ' 

It is well-known that (cf. IB p. 180, Theorem III 4.10.2]) 

Xf'' ^ C([0, t];W''~'''P). (2.3) 
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Let U*''^ be the Banach space of the completion of C^(R'^) with respect to the norm: 

my.,' ■= \\f\L + iiv/ik,p. (2.4) 

For simplicity of notation, we also write 

X*^'" :=Xj''', Y^''' := Yj''' 

and 

W°° := nt,pW^'^, := n^^pY'^'P, X°° := n^^pX*-'', := n^^pU'^'P. 

Let Q. be an open domain of W^. For e No U {oo}, we use = C^(f2) to denote the space of 
all bounded and A:-order continuous differentiable functions with all bounded derivatives up to 
/r-order. For p e (0, 1), let be the completion of <S(R'') with respect to the norm 

ll/lh^/':=ll/llco + l/h^/', 

where || • ||co is the sup-norm and 

I r, l/(^) - /CV)I 

Ify? := sup — — . (2.5) 

x^y \x - y\P 

Notice that C~(R.'') <t By the Sobolev embedding theorem, one has 

W^'PCt^^-f', p>d. 

Let ^ be the class of all linear operators K : W°° W°° satisfying that for each /3 > and 
p> I, 

K : W^ ^' ^ W^-'' is a bounded linear operator, 

and for each /3 e (0, 1), 

mfy<c,,Mfy,, v/e<^^. (2.6) 

A typical example of such an operator is the Riesz transform: 

\y\" 

Indeed, it holds that for any p > I (cf. US), 

l|V/||p - ||(-A)Vllp. (2.7) 

Recalling that for any / 6 SiW'), 

(-A)V(^) = Q r [m-f(x + y)]\y\-'-'dy, (2.8) 
where q > is a universal constant, we have 

(-Ar-ifg) = g(-Ar^f + f{-Ar-g - <^(f,g), (2.9) 

where 



njf := (-A)-^5,/ = lim r fix - y)j:^dy. 



c^(f,g)ix):=c, f {f{x)-f(x + y))(g(x)-g(x + y))\y\-''-'dy. (2.10) 

Jr'' 

From this formula, it is easy to derive that (see lfT6l ). 

Lemma 2.1. Let ^ e SiW') and set ^,(x) := ^{x - z)for z e R'^. Then for any G [1, oo), there 
exists a constant C = C{p, d) > such that for all f e W^'^, 



/ 



-A)HfQ - (-A)V^JI^dz < Cm\ip\\f\\f\\f\\l'l. (2.11) 



For given > 0, / e L°°([0, 1]; W°°) and if e W°°, let us consider the following heat 
equation: 

dfU + Ao(.-A)^u = f, u(0) = (f. 
It is well-known that the unique solution can be represented by 



Jo 



U(t, X) = ri'ifiix) + 



where i'Pf°)t;^o is the Cauchy semigroup associated with Ao(-A)i and given by 



(fiAoy + x) 



(IjP + f2)(./+l)/2 

where c^/ > is a universal constant. By the classical Littlewood-Paley-Stein theory, there exists 
a constant C > only depending on Aq, d,p such that for any / 6 L^([0, 1] x R'') (cf. [[T2l[T6ll ). 



Jo Jo 



V I 'Ptf(s)ds 



ds<C 



Il/(^)ll^d5. 



(2.12) 



We now use the probabilistic technique to extend the above estimate to the more general 
case. Let {L,),^q be a J-dimensional Cauchy process with Levy measure v(djc) = dLxl\xf^^ . It is 
well-known that (cf. L2J) 

^;'V(^) = E^(^ + MLt). 

Let 1? : [0, 1] and A : [0, 1] [0, oo) be two bounded measurable functions. Define 



Tffifix) := Eif 



(2.13) 



By the theory of stochastic differential equation (cf. ^ p.402, Theorem 6.7.4]), one knows that 



d.T^fcp + A{-A)iT,fcp + § ■ ITtfip = 



Now if we define 



then it is easy to see that 



We have 



u{t, x) 



Jo 



f(s, x)ds. 



dtu + A{-Ayu + §-Vu = f, u{0) = 0. 



Theorem 2.2. Let ■& : [0, 1] and A : [0, 1] [Aq, oo) be two bounded measurable 

functions, where Aq > 0. For any p e (l,oo), there exists a constant C depending only on 
Aq, d, p such that for all f e LP([0, 1] x R"'), 

p 



Jo Jo 



fms 



/■ 

Jo 



ds<C WfisWds. 



Proof. Let (L\'-')t^o, i = 1, 2 be two independent copies of Cauchy process {L,)t^Q. By the theory 
of stochastic differential equation (cf. [|2iil5J), one can write 



T'fip{.x) = 



where 'pf"(p(x) := E(p(x + AqL)^') is the semigroup associated with /lo(-A)2, and 

Xr := [ &(r)dr- [ (Air) - Ao)dL^P . 
Jo Jo 



|x - J mdr + J (Air) - Ao)dLi'^ + Ao{Lf^ - Lf ))j = E^f«,^ (x - X, + X,) , 



(2.14) 



Define ^ 

u{t,x):= r r^lJ{s,x + X,)ds. 
Jo 

By (|2.14l) one has 

r r/ZVC'^, ^)d^ = Em {t, X - Xt) 
Jo 

Hence, by Holder's inequality and Fubini's theorem, 



V T;f/(^,x)d5 d?= ||EViia,--x,)|rdf<E ||Vi<a,--x,)|rdr 

Jo p Jo 'Jo 

= E r |KOII^dr = E r V r rtf{s,- + X,)As 
Jo Jo Jo 



< CI 



r'ii/(^.- 

Jo 



+ d5 = C WfisWAs 



Jo 



The proof is finished. □ 
Below we prove a maximum principle for the fully nonlinear equation (11.11) . 

Theorem 2.3. (Maximum principle) Let F(t, x,w,q) : [0, 1] xR^'xR'^ xR — > R be a measurable 
function. Assume that for any R > and all (t, x, w, (?) e [0, 1] x R'' x R'' x R with \w\, \q\ < R, 



< dgFit, X, w, q) < a^j, \^„F\{t, x, w, q) < aR^, 
where aR^i,aR^2 > 0. Let u e C([0, 1]; Cl(R^)) satisfy 

u{t, x) = m(0, x) + I F(s, X, Vu{s, x), -(-A)'m(5, x))ds. 
Jo 

IfFis, X, 0, 0) < 0, then for all t e [0, 1], 

sup u(t, x) < sup m(0, x). 

xeW xeR'' 



(2.15) 



In particular, 

IKOIloo < IKO)|U + I ||F(5,-,0,0)|Ud5. 
Proof. First of all, we assume 



Jo 



F(s, x,0,0)<6< 0. 
Suppose that (12.161) does not hold, then there exists a time to 6 (0, 1] such that 



(2.16) 

(2.17) 
(2.18) 



sup u(tQ, x) = sup sup u(t, x). 

te[0,l] xeR'' 



Let x„ 6 R'' be such that 



lim u(to, Xn) = sup u(to, x) > u(t, x), V(?, x) e [0, 1] x 



We have for any s e (0, to), 

1 



lim u(tQ, Xn) - lim u(tQ - s, Xn) 



< - 

s 

< - \im{u(tQ, Xn) - u(to - S, Xn)) 



(2.19) 



F(S, Xn, Vuis, Xn), -(-A)2m(5, Xn))ds, 
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(2.20) 



and for any h eW', 

< lim - (lim u(to, Xn) - lim m(?o, jc„ - sh)\ 



< lim - lim(M(?o, x„) - u{tQ, x„ - sh)) 



»1 



= lim lim I /z ■ VM(fo, a;„ - ssh)ds 

n->oo Jo 

= lim(/z • Vm(?o, -«n)). 

n— >oo 

In particular, by the arbitrariness of h, we get 

\imVu(to,Xn) = 0. (2.21) 

On the other hand, since for any y eW^, 

u(to, Xn +y)- u(to, Xn) < SUp u(tQ, x) - u(tQ, Xn) ^ 0, 

by (12.81 ) we have 

liin -(-A)5M(fo, Xn) < Q ( \\m[uito, Xn + y) - u(to, Xn)]\yr^~^dy < 0. (2.22) 

Moreover, since by m e C([0, 1]; Cl{W')), 

lim \\V{u{s) - m(?o))IU = 

s-ito 

and 

lim ||(-A)5(m(5) - M(fo))lloo = 0, 

we have by (12.151) . 

1 r'" 1 I 

lim - \\F{s, Vu{s), -{-^)-^u{s)) - F(s, Vuih), -(-A)5M(ro))|Ud5 = 0. 

Hence, by (IZIOl) . (IZITl) and (I2l8]) . 

1 r'" 1 
< lim - lim I F(5, jc„, Vm(?o, Xn), -(-A) 2 1/(^0, x,i))d5 

= lim - lim I F(5, jic„, 0, -(-A)2m(?o, -<^,i))d5 

^^'^ ^ n-^CO Jto-£ 

1 r'" 1 

< lim -lim I [F(5, a:„,0, -(-A)2M(ro, -«„)) - i^(i', 0, 0)]d5 + 5 

= lim lim [a„,£ • ( - (-A)^-u(to, Xn))] + 6, (2.23) 



where 



1 p 1 

an,E-=- I (5^F)(^,x,„0, -r(-A)2M(?o,.x:„))drd5. 



£ 

Leti? := ||(-A)5m(?o)||oo. Noticing that 

< a„,e < aR,i, 

by (12.221 ), (12.231 ) and 5 < 0, we obtain a contradiction. 
We now drop assumption (12.181 ). For 6 <0, set 

us(t, x) = u(t, x) + St. 
1 



Then 

[(5 + F{s, X, Vusis, x), -{-A^usis, x))^ds. 
So, for all re [0,1], 

sup u(t, x) < sup us{t, x) - 6t sup u(0, x) - 5t. 

Letting 5 t 0, we conclude the proof of (|2.16l) . 
As for (12.171) . by considering 

u(t,x) = u(t,x) - I -,0,0)1 |ood5 

Jo 

and using (12.161 ) for u{t, x) and -u{t, x) respectively, we immediately obtain (12.171) . □ 

Next we recall Silvestre's Holder estimate about the linear advection fractional-diffusion 
equation. The following result is taken from [16, Corollary 6.2]. Although the proofs given 
in [[TT| and [fT6l are only for constant diffusion coefficient a(t, x), by slight modifications, they 
are also adapted to the general bounded measurable function a{t, x). 

Theorem 2.4. (Silvestre's Holder estimate) Let a : [0,1] xW' ^ R and b : [0,1] xW' ^ 
be two bounded measurable functions. Let u e C([0, 1]; C^(R'')) satisfy 

u(t,x) = u(0,x) - I {a{-A)^u)(s,x)ds + I {b ■ Vu)(s,x)ds + I f(s,x)ds. 
Jo Jo Jo 

Ifa(t, x) > aa > 0, then for any y e (0, 1), there exist e (0, 1) and C > depending only on 

d, ao, y and HalU, II^IU such that 

sup \u{i)\^, < C(||m|U + ll/IU + l"(0)|^r), (2.24) 

fe[0,l] 

where | - 1<^/3 is defined by ( 12.51) . 

3. Linear nonlocal parabolic equation 
In this section, we consider the following linear scalar nonlocal equation: 

5,M + a(-A)5M + Z7- Vw = /, m(0) = (3.1) 

where a : [0, 1] x R'^ ^ R and : [0, 1] x R'' — > R'^ are two bounded measurable functions. 

An increasing function a; : R'*' — > R'^ is called a modulus function if lim^o (j^{s) = 0. We 
make the following assumptions on a and b: 

(H^'^) Let k £ No, a,b e L°°([0, 1]; C^), and there are two modulus functions coa and ojh such 
that for all t e [0, 1] and x,y e W', 

\a(t,x)-a(t,y)\<cOa(.\x-y\), \bit, x) - b(t,y)\ < cobQx - y\). (3.2) 

Moreover, for some ao, ai > and all (t, x) e [0, 1] x R'', 

ao < a{t, x) < ai. 

We first prove the following important apriori estimate. 

Lemma 3.1. For given p e (1, oo) and k eM, let f e y'''^''' and u e X''''' satisfy that for almost 
all it,x) e [0, 1] xR^, 

d,u(t, x) + a(t, x)(-A)^-u(t, x) + b(t, x) - Vu(t, x) = f(t, x). (3.3) 

Then under (H^'fj), there exists a constant Ck,p > such that for all t e [0, 1], 

WuW^k.p < {\HO%_i^p + ||/|Iy*-i.p j , (3.4) 
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Us(t, x) = Us(0, x) + 



where Ci ^ depends only on uq, ai, \\b\\oo, d,p and cOa, <^b- particular, equation (\3.3\l admits at 
most one solution in X'^'P. 

Proof. Let (pe)ee(o,i) be a family of moUifiers in R'', i.e., Pe(x) = s'''p(s~^x), where p e C^iW^) 
is nonnegative and has support in Bi and Jp = 1. Define 

Usit) := u(t) *Pe, ae{t) := a{t) * p^, b^it) := b{t) * p^, fe{t) := f{t)*Pe- 
Taking convolutions for both sides of (13.31) . we have 

dtu^it, x) + as(t, x)i-Ayueit, x) + bs(t, x) ■ Vu^it, x) = h^it, x), (3.5) 

where 

hit,x) := + aEit,x)i-A)^-Usit,x) - [(a(0(-A)5M(0) * Ps](x) 

+ bs(t, x) • Vus(t, x) - [(b(t) ■ Vuit)) * Pe](x). 

By (|X2l) . it is easy to see that for all s e (0, 1) and t e [0, 1] and x,y e W', 

\as(t,x) - ae(t,y)\ ^ C0a{\x-y\), \bs(t,x) - bs(t,y)\ ^ cotQx - y\), (3.6) 

and 

\ae(t, x) - a(t, x)\ < coais), Ibsit, x) - b(t, x)\ < C0b{s). 
Moreover, by the property of convolutions, we also have 

lim r \\hM-m\\';dt = 0. 
Jo 

Below, we use the method of freezing the coefficients to prove that for all t e [0, 1], 

Wusiy^.p < C (\\us(0)\\i_i p + \\hs\U^o.„y (3.7) 

where the constant C is independent of s. After proving this estimate, (13.41) with k = I immedi- 
ately follows by taking limits for (13.71) . 

For simplicity of notation, we drop the subscript s below. Fix 6 > being small enough, 
whose value will be determined below. Let ^ be a smooth function with support in Bg and 

UWp = 1. ForzeR'', set 

ax):=ax-z), i^=a(f,z), ^'^it):=bit,z). 
Multiplying both sides of (13.51) by 4^ we have 

dMz) + A"^i-A)ku^z) + ^\ ■ V(m4) = si, 

where 

:= (i^ - a)(-Ar^u^, + A%-Ar-iu^,) - (-Ay^uQ 



-AIM 



+ {&l-b)-V{uQ + ub-V^, + H,. 
Let T^""' be defined by (I2l3l) . Then w^, can be written as 

uUt,x) = r^l''ku{0)Q{x)+ r rfY'"^g[{s,x)As, 

Jo 

and so that for any T e [0, 1], 

r mu^mrpdt<2P- 

Jo 



( r-T 



IIVT;J''~-(M(0)^,)ll^d? + 



V I T;^'^gl{s)ds 



dt 
p ) 
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For h{T, z), by (IXT?1) and (O), we have 

\\VT'j;''ku{0)Q\\ld.t = \\V'P^\u{0)Q\\Ut 



T 



dt 

p 



<C||KO)^,|i;_, (3.8) 

where the last step is due to [fT4l p. 96 Theorem 1.14.5] and (12.21) . Thus, by definition (12.11) . it is 
easy to see that 

r hiT,z)dz <c f wumzW' . dz < c(iim(0)ii; . UK + iim(0)iip^ii; , J. 

For hiT, z), by Theorem l2.2[ we have 

UisWpds < C \m't - a){{-Ayums)\\';ds 
Jo 

+ C r |U;'((-A)^i<^,-(-A)^(i<4))(5)||^d5 
Jo 

+ c r ii((i^,^-Z7)-v(i.4))(^)ii;;d^ 

Jo 

+c r \\{ub-vQ{s)\%ds+c r ii/ia^)ii^d^ 

Jo Jo 
=: /2i(r, z) + /22(r, Z) + /23(r, Z) + /24(r, Z) + /25(r, Z). 

For /2i(r,z), by (IJ!6l) and ||^||p = 1, we have 

r hi{T,z)dz<CcoP{6) r r ||((-A)^M4)(^)||j^d^dz 
Jr"* Jw' Jo 

\\{-Ayu{s)%ds < Ct^^(5) ||VM(5)||;d5. 
Jo 

For 722(7", z), by (12.1 ID and Young's inequality, we have 

r 722(r,z)dz<Cai r r ||((-A)^M4-(-A)kM^J)(^)ll^dzd5 

Jr"* Jo Jr"* 

\\u{s)\%ds + c{ \\u{s)\\f\\Vu{s)\\fds 
Jo Jo 

< C J^^ MsWpAs + WVuisW^ds. 

For 723 (r, z), as above we have 

j^^ 723(r, z)dz < Ccolid) I \\Vu(s)\\Pds + IWa'p Il"(*)li;;d5j . 

Moreover, it is easy to see that 

r 724(r,z)dz<C||Z;||^||V^Ii;^ r IK^)||^d5, 
Jr'' Jo 



r 725(r,z)dz<c r ii/i(5)iii^d5. 

Jr"* Jo 
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Combining the above calculations, we get 



f \\Vu(s)\\Pds= f f \\Vu(s)-a';dzds 
Jo Jo Jr"* 

<2P-' [ [ \\V(uQ(s)\\Pdzds + 2P-'\\Va';, [ MsW^ds 
Jo Jr'' Jo 

< C\\um\l, ^ + (\ + C(coP{6) + co^m) J misWpds 

+ C [ Ms)\\Pds + C [ WKsWpds. 
Jo Jo 



Choosing 60 > being small enough so that 

1 
4 

we obtain that for all T 6 [0, 1], 



CicoPi6o) + coli6o))<-, 



[ \\VuisWpds<C\\u(0)r +C f MsWpds + C f ms)\\';ds. (3.9) 
Jo p'^ Jo Jo 

On the other hand, by (13.51 ), it is easy to see that for all t e [0, 1], 



< ci|m(0)||^ + c iiVM(^)ii;^d^ + c msW^ds, 

Jo Jo 
which together with (13.91) and Gronwall's inequality yields that for all t G [0, 1], 

sup|K^)||^+ r WVuisXds^ciwumi' ^ + f WHsWpds]. (3.10) 

i€[0,f] Jo \ i''P Jo I 

From equation (|3.3I) . by (12.71) we also have 

f\\dMsXds<c(\\a\t f \\{-Ar^u{sXds + \\b\t f WVuisX^s + f \\h{s)\\';,ds] 
Jo \ Jo Jo Jo / 

< C |(||a||^ + ||Z7|P:) WVuisXds + WHsX^s^ , 

which together with (13.101) gives (13.71) . and therefore (|3.4I) with k = I. 
Let us now estimate the higher order derivatives. For « = 1, 2, ■ ■ ■ , A:, let 

w^"\t,x) := V"u{t,x). 

By the chain rule, we have 

dtw^"^ + a(-A)5w(") + b ■ Vw^"^ = h^"\ 

where 

n I 

Thus, by (|3.4I) with k= \ and the assumptions, we have 

' n 

||V"m(0)||i_i^^ + W^"'u\\., + ||V"/IIyO,. 



< c 
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which implies that 

WuWj^n+i.p < C^||M(0)||„+i_i p + IImIIy';" + ll/lly;"-")- 
By induction method, one obtains (13.71) . □ 
Now we prove the existence of solutions to equation (13.11 ). 

Theorem 3.2. Let p > \ and e N. Under for any ip e W^"^''' and f e Y'^-^^p, there 

exists a unique u e X*^'^ with m(0) = (p solving equation ( 13. 7 P . 

Proof. We use the continuity method. For /I e [0, 1], define an operator 

Ua := dt + ia(-A)5 + . V + (1 - /l)(-A)5. 

By (12.71 ). it is easy to see that 

Ua-.^^''' ^Y^-^'P. (3.11) 

For given (p e W'^''/''^, let X^''' be the space of all functions u e X'^'P with m(0) = (p. It is clear 
that X^''' is a complete metric space with respect to the metric || • ||x*.p. For A = and / e Y'^'^'P, 
it is well-known that there is a unique u e X^'^ such that 

Uou = d,u + (-A)^a = /. 
In fact, by Duhamel's formula, the unique solution can be represented by 

u{t,x) = r]<p{x)+ { 'PlJis,x)ds. 
Jo 

Suppose now that for some Aq 6 [0, 1), and for any / 6 Y*""''^, the equation 

Uao^^ = f 

admits a unique solution u e X^''. Then, for fixed / e Y*"^^'^ and A e [Aq, 1], and for any 
u e X^'^, by (|3.11l) . the equation 

UA,w = f + (UA,-UA)u (3.12) 

admits a unique solution w e X^''. Introduce an operator 

w = QaU. 

We want to use Lemma [3TT] to show that there exists an £ > independent of Aq such that for 
allie [Aq,Aq + £], 

Qa : X'/ ^ X'/ 

is a contraction operator. 

Let Ml, M2 s X^'^ and w; = Qau^, i = 1, 2. By equation (|3.12l) . we have 

Ua,(wi - W2) = (Ua, - Ua){u, - u^) = (^ - A)(ia - 1)(-A)5 + b ■ V)(mi - U2). 
By (l3ll) and (12771) . one sees that 

WQaUi - QAUiiy-P < CkJAo - A\ ■ ||((a - 1)(-A)5 + b ■ V)(mi - M2)||y*->./' 
< CoUo - ^1 ■ 11^1 - "iIIy*." < Cq\Ao - A\ ■ \\Ui - M2||x*.P, 
where Co is independent of A, Aq and U\,U2- Taking £ = 1 /(2Co), one sees that 

Qa : X^P ^ X^P 

is a 1 /2-contraction operator. By the fixed point theorem, for each A e {Aq, Aq + e], there exists 
a unique u e X^'' such that 

QaU = M, 
12 



which means that 

Uau = /. 

Now starting from A = 0, after repeating the above construction [-] + 1 -steps, one obtains that 

for any / e Y*^-'-'', 

U,u = f 

admits a unique solution u e X^''. □ 

4. Quasi-linear nonlocal parabolic system 
Consider the following quasi-linear nonlocal parabolic system: 

dfU + a(u,'RaU)(-A)^u + b(u,'RbU) ■ Vu = f(u,'RfU), (4.1) 
where u = (u\ - ■ ■ , a™) and 

a{t, X, u, r) : [0, 1] x x R"" x R*^ ^ R, 
b{t, X, u, r) : [0, 1] x R^ x R"" x R*^ ^ R^, 
fit, X, u, r) : [0, 1] X R^ X R"' x R'^ ^ R"", 
are measurable functions, and 

Here we have used that 'R^u = f^'Ju-', similarly for "RhU and "RfU. 
The main result of this section is: 

Theorem 4.1. Suppose that a(t, x, u, r) > gq > 0, and for any R > 0, 

a,b,fe L'^iiO, 1]; cr(R^ x B"^ x 4)), (4.2) 

a, be L'"([0, 1]; C,J(R^ x B"j} x R'^)), (4.3) 

where B'^ denotes the ball in R™ with radius R, and for each j e No, there exist Cfj, jj > and 
hj e (Li n L'^)(W') such that 

Wifit, X, u, r)\ < Cfj\u\(\u\->'' + 1) + hj(x), (4.4) 

and for some Cf > 0, 

(u, fit, X, u, r))R„, < C/(|Mp + 1). (4.5) 
Then for any (p e W", there exists a unique u e X°° solving equation A4.1\l . Moreover, 

sup Mt)\\i < ^''^(ll^lli + Cf). 

(£[0,1] 

Proof. First of all, for any 'R e ^ and u e X'''^, by the boundedness of R in Sobolev space W'''^, 
one has 

{t,x)^'Ru{t, x)eX'''P. 
Thus, by (14.21) and the chain rules, one sees that for any u e X'", 

(t, x) ^ a(t, X, u{t, x), Raiiit, x)) 6 L°°([0, 1]; C^), 
{t, x) ^ b{t, X, u{t, x), Rbuit, x)) 6 L"([0, 1]; C"), 

and by (1441) . 

{t, x) fit, X, u(t, x), Rfuit, x)) e y". 

Set uo{t,x) = 0. By Theorem |3.2[ we can recursively define w,, e X~ by the following linear 
equation: 

dtUn + a(Un-l,KUn-l)(-AyUn + b{Un-l,RbUn-\) ' Vm„ = f{Un^\,RfU„^l) (4.6) 
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subject to the initial value u„{0) = (p e W". 

We first assume that 70 = 0. For simplicity of notation, we set 



bn{t, X) 
fnit, X) 



- a(t, X, Un-l(t, X), 'RaU„-i(t, x)), 
= b(t, X, U„-i{t, X), KbUn-iit, x)), 

= f(t, X, u„^i{t, x), 'RfU„^i{t, x)). 



By the maximum principle (see Theorem I2.3D and in view of yo = it is easy to see that 

llWnWIU < ||m„(OIIoo + J 11/(5, •, Un-l(s, ■), KfUn-lis, ■))IUd5 

< I|W«(0)|U + ^(C/,oll"«-l(5)IU + I|/^olloo)d5 

< WfWoo + ll^olloo + Cffl J ||M„_i(5)||ood5, 

which yields by Gronwall's inequality that 

sup ||m„(OIIco < e^^"(llvlU + WhoWoo) =: Ko. (4.7) 

tern] 

By Theorem 12 .41 and (|2.6I) . there exist /3 e (0, 1) and constant ^1 > depending on ^0 such that 
for all n eM, 

sup \Un(t)\^^P + sup \KUn(t)\^P + SUp \'RbUn(t)\<^P < Ki. 
te[0,l] fe[0,l] f€[0,l] 

Thus, by (14.31) we have 

\a„{t,x)-an{t,y)\ < ||V.,a||i~ |x - j| + i^i(||V„a||i» + ||V„a|L~ )|x - (4.8) 

\b„(t,x)-bn(t,y)\ < \\VML-Jx-y\ + K,(\\VML-^ + Wub\\L-J\x-yf, (4.9) 

where || • denotes the sup-norm in L'"([0, 1] x R"' x 5™ x R*^), 
For^ = 0, 1,2, • • • , set 

wf\t,x) := V^u„(t,x). 

By the chain rule, we have 

5,wf + a„(-A)5w« + Z7„ • Vwf = , 
where ^ = fn and for > 1 , 

By (1481) . (1491) and Lemma[3TIl we have for all p > 1 and t e [0, 1], 

llwf 11^,,. < Qp(||VVlli_i,p + ll^f 11^0,,), (4.10) 

where Ck^p is independent of n. 

For = 0, by (OOl) . (l44l) and (lO) . we have 

iMoii^+ r iM5)ird5 < ciMi; , +c r ii/„(5)ii^d5 

Jo ' p'^ Jo 

< Cll^li; ,^^^ + (Cf,o(Kl' + l)\\u„.,{s)\\p + \\ho(s%yds 

<C|Mi;, +C r \\Un-ds)\\Pds + C [ \\ho(s)\\Pds. 
p''' Jo Jo 
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By Gronwall's inequality, one gets 

sup sup ||M„(Ollp < Cp, 

nen re [0,1] 

and therefore, for all p > 1, 

SUp||M„||xi.P <Cp. 

neN 

Now for any k = 1 , 2, ■ ■ ■ , since by the chain rules, g^^^ only contains the powers of all derivatives 
up to fc-order of m„, 'RaUn, 'RhUn and RfUn, by induction method and using Holder's inequality, it 
is easy to see that for all e N and p > I, 

sup \\Un\y.i' < Ck.p. (4.11) 

Below we write 

Wn,m(t, X) := Unit, x) - U,„(t, x). 

Then 
where 

By Lemma [3T| again, we have for all ;? > 1 and t e [0, 1], 

l|l^«,mllx'-'' ^ C||^„_mllY°'''- 

Here and below, C > is independent of n, m. Using (14.1 II) and (|4.3I) . we have 

ll,?n,mllY°-'' ~ /kiIIyJ'-'' W'^n ~ '^mllYj'-'' H^" ~ ^mllYj"-'') 

< C(||V„/|L» + ||VJ|L» + ||V„a|L» + ||V,a|k» 

+ I|V„Z7||l~ + \\VM\l'^ )l|w„-l,,n-lllYO-f- 

Hence, 

sup < C I \\w„.i,„^iis)\fds. 

ie[0,f] Jo 

Taking sup-limits and by Fatou's lemma, we obtain 

lim sup ||w„,™(5)|H^ < C I lim sup ||w„,,„(5)||^dr. 

"'"'^"^ ie[0,r] Jo ie[0,r] 

So, 

lim sup ||w„,„,(*)ll^ = 0, 

which together with (14.111) and the interpolation inequality yields that for all ^ e N and p > I, 

lim sup = 0. 

Thus, there exists au e X°° such that for all e N and p > I, 

lim sup \Kis) - uis)\\l = 0. 



«,m— ►oo 



Taking limits for (14.61) . one sees that u solves equation (14.11) . 

Now we want to drop yo = and assume (14.51 ). For R > 0, letXR 6 C^(K.'') be a nonnegative 
cutoff function with xr(u) = 1 for \u\ < i? and;^K(M) = for |m| > + 1. Set 



fR{t, X, u, r) := f{t, X, u, r)xR{u) 
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Let Ur g X~ solve 

dtUR + a{uR,%aUR){-^yuR + b(uR,'RtUR) ■ Vur = fR{uR,'RfUR). 
Noticing that by (|Z9l) . 

- -2 
2((-A)5Ms,Mfi)R" = (-A)2|Mr| + &{Ur,Ur), 

we have 

ldt\uR^ + a(MR,KMi?)(-A)5|M«p + biuR,'RhUR) ■ V\ur\^ 

= 2{uR,fR(uR,'KfUR))^,.,-a(uR,'RaUR)&(uR,UR) < 2Cf(\uR\ +1). 
Thus, by the maximal principle, we have 

WuRimi < ii<^iii + CfJ^ (\\uR(s)\\i + i)ds, 

which implies that for all R > 0, 

sup WuRimi < e'^fiMi + Cf). 

f6[0,l] 

The proof is finished by taking R := + C/)]'^^ □ 

5. Fully nonlinear and nonlocal equation: Proof of Theorem I1.1I 

The following lemma will play a key role in proving the existence. 

Lemma 5.1. Leta e L'"([0, 1];C^(R^)) be bounded below by ao >Oandb e L°°([0, l];Cl(W')). 
Let u : [0,1] xW' ^ belong to X^'P for some p > 1 and satisfy 

dtu = ai-A)~'^nu + b-{Vu- (Vu)'), (5.1) 
where □ : = divV - Vdiv. Then we have 

\\u\y. + \\u\y,, < c(iii/(0)iii,p + iit/(0)lli,p), 

where U ■.= Vu- (Vu)'. 

Proof. By equation (|5.1I) . one sees that 

dtU = -a{-A)^u - a(-A)"5VdivM + b ■ U, 

and 

dtU = -a(-A)5 [/ + b-VU + (Vb) ■ U - [(Vb) ■ U]' + A, 

where 

A := (Va)' ■ (-Ay-^nu - (i-Ay-^nuY • Va. 
By Lemma [Xn there exists a constant C > such that for all t e [0, 1], 

\\u\y.p < C||M(0)||i„i p + C\\a(-Ay^Vdwu\\^o.p + C\\b ■ U\\^o.p 

< C||m(0)||i„ 1 „ + C||a||oo||divM|Lo,P + C\\b\\oo\\U\Up, 

IL^ JL ^ 

and 

\\u\y.p < ciit/(0)lli-i,, + cKVb) -u + u- (Vby+A\\^o,, 

< C||f/(0)lli-i p + C(||Vfl|U + ||VZ7|U)||Vm||yO,.. 
In particular, for all t e [0, 1], 

IKOII^ + J Ms)\\lds<C\\um\l^ + C J WdiYuisXds + C J \\U(sXds 
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< C\\u(0)f + Ctl sup \\divuis)\\P + sup \\U(s)\\A , (5.2) 

and 

mtWp + \\u{s)\\'; ds < cwumiip + cj^ Ms)\\[ds. (s.s) 

On the other hand, noticing that 

divDw = 0, 

we have 

dfdi\u = Va ■ {-Ay^nu + diy(b • U). 

Hence, 



||divM(OII,, < l|divM(0)|L + llValU \\(,-Ar^Du(s)\Lds 



a\\oc I 
Jo 



+ \\b\\^ \\VU(s)\\pds + \\Vb\\^ ms)\\pds 
Jo Jo 

< C(||divKO)||p + ||Wlli,p) + Cn . (5.4) 

Now substituting (|53] ) and (|54l ) into (|?!2l) , we obtain that for all ? G [0, 1], 

Mm';, + iii<(^)iiy^ < Co(\\um\ip + \\um\ip) + c^tj^ 11^^(^)11^',/^, 

where Co, Ci are independent of ||m(0)||i,p and ||i7(0)||i_p. Choosing to := l/(2Ci), we arrive at 
sup lKOIi;^+ ^11^^)11^/^ <2Co(lKO)ir + ||[/(0)lll\,). 

te[0,to] Jo 

So, for some C2 > 0, 

\\u\\' + \\u\f < C2(\\um\';^ + \\um\lp). 

In particular, 

f (ms)\\', + \\u{s)\f )ds < C2(\\umf + \\um\'; ). 

Thus, there is at least one point sq e [2fo/3, to] such that 

Mso)\\ip + \\uiso)\\ip < ^(ii"(0)iiL + mo)\\Q. 

Now starting from as above, one can prove that for the same to, 

3(^2 

M- + So)\\' + \\U(- + So)\\' < C2(\\u(so)\\lp + \\Uiso)\\l^) < -^(ll"(0)||[^^ + WUmi'Q. 

'0 '0 '0 

Repeating the above proof, we obtain the desired estimate. □ 
We are now in a position to give 

Proof of Theorem 1 We divide the proof into three steps. 

(Step 1). In this step we consider the following fully non-linear and nonlocal parabolic 
equation: 

dfU = F(t, X, Vu, -(-A)'m), u{0) = (f. 
As introduced in the introduction, let 

= (-A)~^divw. (5.5) 
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For any ip e U°° = C\k^p"U^'P, where U'' '' is defined by (12.41) . by Theorem |4.1[ there exists a unique 
w e X°° solving the following parabolic system: 

d,w = -idgF){w,'Rw){-A)^w + {V„F){w,nw)Vw + V^F{w,nw) 
subject to w(0) = V^. Define 

u{t, x) := ip{x) + I F{s, X, w(s, x), Kwis, x))ds 
Jo 

and 

h(t, x) := Vu(t, x) - w(t, x). 

Then we have 

dth = (dgF)(w,'Rw)(V'Rw + (-A)2w) + (V^F)(w,'Rw)((Vwy - Vw) 
= (d^F)(w,'Rw)(-Ar-^nh + (V„F){w,'Rw)(Vh - (Vh)') 
subject to h(0) = 0, where □ := divV - Vdiv. By Lemma [STTl we have 

/j = => w = Vm. 

Thus, by (153]) . 

dtu(t, x) = Fit, X, Vu(t, x), KVuit, x)) = Fit, x, Vuit, x), -i-A^uit, x)). 
By the maximum principle (see Theorem [23]), we have 



Jo 



IKOIU<IMU+ I ||F(5,-,0,0)|Ud5. (5.6) 

In particular, u e C([0, 1]; U"). 

(Step 2). Now we consider the general case. Set uoit,x) = 0. Let m„ e C([0, 1];U'") be 
defined recursively by the following equation: 

dtUn = Fit,X,Un-\,VUn,-{-AT-Un), Un{G) = (fi. (5.7) 

By (|5^ and (fL6l) . we have 



Jo 



"«(Olloo < \\(f\\oo + \\Fis, •, Un-lis, ■), 0, 0)||ood^ 



< ll^lloo + Ko i\\Un-lis)\\c^ + l)ds. 

Jo 

By Gronwall's inequality, we get 

WuML < e^'iMoo + Ko) =: Ko. (5.8) 
On the other hand, by taking gradients with respect to x for equation (15.71) . we have 
dtVun = -dqFit, X, u„-\,Vun, -(- A) 2 M„)(-A) 2 Vm„ 

+ V^Fit, X, U„-i,Vu„, -i-A)^Un)V^Un 
+ duFit, X, Un-i,VUn, -(-A)5m„)Vm„_i 

+ V^Fit, X, Va,j,-(-A)5M„). 
By the maximum principle again and (11.91 ). (11.101) with y^^ i = 0, we have 

I|Vm„(OIIoo < l|V(^|U + I \\duFis,X,Un-l,'VUn,-i-A)^Un)'VUn-l\\oods 

Jo 

+ ( \\V_^Fis,x,u„-i,Vu„,-i-A)^Un)\\oods 
Jo 
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< IIV^IU + C r (\\VUn-lis)\\oo + \\VUnis)\\oo + l)ds. 



where C is independent of n. By Gronwall's inequality, we get 

sup sup ||Vm„(OIIoo < +00. (5.9) 

n te[0,\] 

Moreover, by (fLSl) . (fLQl ). (11.101) . (15^ . Theorem|23]and Lemma[2?IJ as in the proof of Theorem 
14. 1[ we have for all p > 1, 

||Vm„||^i,p < C^||V^||i_ip + ||Vm„_i||yO,p + ||Vm„||yO,,. + 

which implies by Gronwall's inequality that 

sup ||Vm„||xi.p < +00, (5.10) 

n 

and furthermore, for all ^ e N and p > I, 

sup ||Vm„||x*.p < +00. (5.11) 

n 

This together with (15.81 ) gives 

sup sup ||m„(OIIu*j' < +<^- (5.12) 

n re[0,l] 

(Step 3). Next we want to show that m„ converges to some u in C([0, 1]; tf '^). For n, m e N, 
set 

Vn,mi,^i '•— Unit, X) — Ujnit, X). 

Then 
where 



id^F){Un-i,Vun, -(-A)5(*v„,^ + uj)ds, 



Jo 





1 



{V^F){u„-U^{sVn,m + Mm), -(-A)5M„)d5, 



1 

fn,m-= I (5„F)(5V„_i,^_i + M^_l, VM,„,-(-A)5Mjd5. 

Jo 

By the maximum principle, we have 

l|v„,m(OIU < C I ||v„-i,m-i(^)|Ud^, 

Jo 

and by Gronwall's inequality, 

lim sup ||v„,„,(OIU = 0. (5.13) 

"■'"-^«'fe[0,l] 

On the other hand, by Lemma [ZTI and (15.121) . we may derive that for all t e [0, 1], 
and so, 

lim ||v„,,„||xi.P = 0. 

n,m— >c» 

This together with (|5.11l) . the interpolation inequality and (|5.13l) yields that for all A; e N and 

p>l, 

lim sup ||v„,m(Ollu*p = 0. 

«.'«-^'»f€[0,l] 
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Thus, there is a m e C([0, 1]; U°°) such that for alU e N and p > I, 

lim sup \\Un(t) - m(OIIu*;' = 0- 

The proof is complete by taking limits for approximation equation (15.71) . 
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